Series expansion studies of random sequential adsorption with diffusional relaxation 
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We obtain long series (28 terms or more) for the coverage (occupation fraction) 9, in powers of 
time t for two models of random sequential adsorption with diffusional relaxation using an efficient 
algorithm developed by the authors. Three different kinds of analyses of the series are performed 
for a wide range of 7, the rate of diffusion of the adsorbed particles, to investigate the power law 
approach of 9 at large times. We find that the primitive series expansions in time t for 9 capture 
rich short and intermediate time kinetics of the systems very well. However, we see that the series 
are still not long enough to extract the kinetics at large times for general 7. We have performed 
extensive computer simulations employing an efficient event-driven algorithm to confirm the t ' 
saturation approach of 9 at large times for both models, as well as to investigate the short and 
O l_ intermediate time behaviors of the systems. 
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INTRODUCTION 

Random sequential adsorption (RSA) |IJ is an irreversible process which particles are deposited randomly and 
consecutively on a surface. The depositing particles, represented by hard-core extended objects, satisfy the excluded 
volume condition where they are not allowed to overlap. The exclusion of certain regions for further deposition 
attempts due to the adsorbed particles leads to a dominant infinite-memory correlation effect where the system 
approaches partially covered, fully blocked stage at large times. However, this picture is altered when the diffusional 
relaxation is introduced ■ Privman and Nielaba || have shown that the effect of added diffusional relaxation in 
the deposition of dimer on a ID lattice substrate is to allow the full, saturation coverage via a ~ t^ 1 / 2 power law at 
large times, preceded by a mean- field crossover regime with the intermediate ~ t^ 1 behavior for fast diffusion. 

Series expansion is one of the powerful analytical methods in the RSA studies [^-[To[ . Long series in powers of time 

St have been obtained, reminiscent of series expansions in equilibrium statistical mechanics, by using a computer 
Recently, the authors [[To[ have proposed an efficient algorithm for generating long series for the coverage 9 in powers 
■ of time t based on the hierarchical rate equations. 

The present work is to study the time-dependent quantity 9 for one-dimensional models of RSA with diffusional 
relaxation, both analytically and numerically. It will be seen that even though relatively long series have been 
^ ' obtained, we are still unable to extract the kinetics of the systems at large times for general 7 due to long, rich 
^ . transient crossover regime that the series must describe. Extensive computer simulations are performed to confirm 
the i -1 / 2 power law approach of 9, where we have employed an efficient event-driven algorithm. The remainder of 
this paper is organized as follows. Section | introduces two related models. Details of series expansion are explained 
in Section Analyses of the series can be found in Section [II . Monte Carlo results are presented in Section [V and 
finally Section ^ contains the summary and conclusions. 



I. THE MODELS 



Two models have been studied in this work. We start with an initially empty, infinite linear lattice. Dimers are 
dropped randomly and sequentially at a rate of k per lattice site per unit time, onto the lattice. Hereafter we set k 
equal to unity without loss of generality. If the chosen two neighbor sites are unoccupied, the dimer is adsorbed on 
the lattice. If one of the chosen sites is occupied, the adsorption attempt is rejected. One of the simplest possibilities 
of diffusional relaxation in this dimer adsorption process is that the adsorbed dimer is permitted to hop either to 
left or right by one lattice constant at a diffusion rate 7 from the original dimer position, provided that the diffusion 
attempt does not violate the excluded volume condition. This model has been initiated and studied by Privman and 
Nielaba 0]. We refer this model as the dimer RSA with dimer diffusion or diffusive dimer model. 

A second possibility is that an adsorbed dimer is allowed to dissociate into two independent monomers; each 
monomer can diffuse to one of its nearest neighbor sites with a diffusion rate 7, provided that the diffusion attempt 
does not violate the excluded volume condition. This model bears a strong resemblance to the former model and 
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is exactly solvable when 7 = 1/2 |1J]. We refer this model as the dimer RSA with monomer diffusion or diffusive 
monomer model. Interestingly enough, the special case of the diffusive monomer problem with 7=1/2 can be mapped 
to the diffusion-limited process 

A + A^ inert, (1) 

which is known as one-species annihilation process ]l3| . This model has been solved exactly by a number of researchers 
Km. We observe that when 7 = 1/2, the effect of a dimer deposition attempt in the diffusive monomer model 
corresponds to two diffusion attempts of A in an adjacent pair of A of the .A + .4 — > inert process. The time-dependent 
quantity coverage 9(t) (fraction of occupied sites) for the diffusive monomer model with 7 = 1/2 is given by 

9{t) = 1 - exp(-2t)/ (2i) ! (2) 

where I n (z) is the modified Bessel function of integer order n. 

II. SERIES EXPANSIONS 

To illustrate how series expansions are performed, we note that the first few rate equations for the dimer and 
monomer diffusive models are 



and 
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respectively, where P(C) denotes the probability of finding a configuration C of sites specified empty 'o' or filled 
Unspecified sites can be occupied or empty. Here we have taken into account the symmetries of a configuration under 
all lattice group operations. For the one-dimensional configurations, we just need to consider the reflection operation 
only. 

Let C denote a particular configuration of interest, and Pc a = P(C Q ) the associated configuration probability. Pc a 
is expected to be a well behaved function of time t, so one can obtain the Taylor series expansion with the expansion 



point at t = 0, Pc„{t) = ]C Pc € } n H n /nl, with the nth derivative of Pc a given by 

71=0 

(n) d n P Co {t) 
PcV - 



dt" 



(10) 



Let Gi denote the set of new configurations generated in the calculation of the ith derivative of Pc , and G\ the 
corresponding jth derivatives of the set of configurations. We observe that Gg _1 1 G" -2 , . . ., G^ l _ 1 (determined at 
the (n — l)th derivative), Gq~ 2 , G™~ 3 , . . ., G°_ 2 (determined at the (n — 2)th derivative), . . ., G° are predetermined 
before calculating the nth derivative of Pc a ■ In the calculation of nth derivative of Pc a , we determine systematically 
Gq, G™~ , . . ., G\_ 1 , G°, by recursive use of rate equations. This algorithm is efficient since each value in G"~ l , 
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< i < n and the rate equation for a configuration C is generated once only. However, this algorithm consumes the 
memory quickly as a result of storage of intermediate results. 

The computation of the expansion coefficients makes use of the isomorphism between a lattice configuration and 
its binary representation if we map an occupied (empty) site to 1 (0). The data structures used to represent Eq. (||) 
and Eq. (Q) are depicted in Fig. |l|. A node for a configuration C is characterized by its four components; (i) the 
representation of C in the computer, (ii) a pointer to the derivatives of Pc, Pc for n — 1,2,3,..., (hi) the highest 
order of derivative h of Pc obtained so far, and (iv) a pointer to a linked list of nodes of configurations ('children') 
which appear in the right hand side of the rate equation for Pc- The linked list contains the associated coefficients for 
each 'child'. The variable h is used so that we know the values of Pc^ where 1 < n < h have already been calculated 
and can be retrieved when needed. All pointers to the configuration nodes generated during the enumeration process 
are stored in a hash table or a binary tree to allow efficient checking of the existence of any configuration. Use of the 
algorithm and data structures allows us to obtain coefficients up to t 31 and t 27 (presented in Appendix [A|) for P(o, t) 
of the diffusive dimer and monomer models, respectively. 



III. ANALYSES OF SERIES 

Analytically, we are interested in confirming the power law approach of t~ 1 / 2 of the coverage 9 at large times for 
both diffusive dimer and monomer models through the unbiased and biased analyses of the series. The unbiased 
analysis does not fix the saturation coverage of the system, while the biased analysis assumes the saturation coverage 
to be the value 1. For the unbiased analysis, let 9(t) = 1 — P(o, t) be the time dependent coverage. Let assume that 
at very large times t the coverage 9 satisfies the equation 

9(t) = 9 c -^, (11) 

where 9 C is the saturation coverage. A(t) is assumed to be a function of t, which tends to a constant value as t — > oo, 
and 6 is the exponent that characterizes the saturation approach (we expect to obtain S — 1/2 from the analysis of 
the series for all 7 values). Writing t = A(t) 1 / 5 (6' c — 9)~ 1 / s , we see that if we perform a DLog Pade Jl7|-Jl9[ analysis 
to the inverted series t = t(9), where 

'h^.i.-L., (i 2) 

then the power law of Eq. ([ll]) implies a simple, isolated pole of 9 C with an associated residue of -1/5. Fig. I shows 
the plot of the inverted series t versus 9 for the 28-term series with 7 = 1/2 for the monomer diffusive model. 

For the diffusive dimer problem, the closest real pole to the value 1 (the expected saturation coverage) for [16,15], 
[15,16], [15,15], [16,14], and [14,16] Pade approximants are shown in Fig. [| with the corresponding saturation ex- 
ponents S displayed in Fig. [|. Similarly we form [14,13], [13,14], [13,13], [14,12], [12,14] Pade approximants for the 
diffusive monomer problem, where the results are displayed in Fig. |^ and Fig. ^. Comparing the graphs for these two 
models, the diffusive dimer series give a better convergence of 9 C and S against 7 than that for the diffusive monomer 
series generally, presumably due to the fact that the coefficients of the series of P(°, t) alternate in signs in the former 
model. For small 7 values (7 < 5), the estimates for 9 C and <5 are unstable — different Pade approximants do not 
agree with one another. The series with 7 = describes a pure lattice RSA behavior where the system approaches 
the jamming coverage exponentially. Hence we expect the confirmation for power law of Eq. ( pd| ) is interfered by the 
exponential behavior of the series when 7 is small. For 7 > 5, there are physically favorable estimates for 9 C and 5 
where 9 C = 1.00 ± 0.05 and 5 = 1.0 ± 0.1 for 10 < 7 < 20, for both models. These results are the manifestations of 
the transient regime of i _1 approach to saturation. 

The distribution plot of the poles and zeros in the vicinity of (1,0) is displayed in Fig. [7] for the [14,13] Pade 
approximant for the 28-term series with 7 = 1/2 for the diffusive monomer model. We see that the real pole closest 
to (1,0) is not distinguished and isolated from the nearby poles and zeros. This explains the difficulty of unbiased 
analysis that the intermediate crossover effect masks the power law approach at late stages. 

We also perform biased analyses for the series. This series analysis have been used by Jensen and Dickman |2^] to 
extract critical exponents from scries in powers of time t. We define the ^-transform of f(t) by 

F[f(t)}=t±kif. (13) 

If / ~ At~ a for some constant A, then F(t) — > a as t — + 00. We consider the exponential transformation 
1 - e- bt 
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which proved to be very useful in the analysis of RSA series [|6|,[10 23 1. This transformation involves a parameter b 
which cannot be fixed a priori is then followed by the construction of various orders of Pade approximants to the 
z-series. Crossing region is then searched for in the graphs of a versus b, the transformation parameter. 

To illustrate this biased analysis, we take the saturation coverage 9 C to be 1 and choose f(t) to be P(°, t) = 8 c — 6(t). 
Since we expect P(o,t) ~ r 1 / 2 for large times t, specifically we have formed [14,13], [13,14], [13,13], [14,12], [12,14] 
Pade approximants to the z-series for the 28-term series with 7 = 1/2 for the diffusive monomer model. We find that 
the estimates for 5 is 0.5061(5), for 0.45 < b < 0.50, as we can see from Fig. |[ Thus the exact analytical function 
of exp(— 2i)Io(2t) serves as a useful guide of this analysis, where the exponent deviates from the value 1/2 by only 
about 1%. 

Given a value of 7, we obtain the corresponding estimates of S from the first convergence of all Pade approximants 
by locating the crossing region. The results of 5 estimates for several values of 7 are presented in Table |. The 
corresponding uncertainties for 5 which reflect the variation of 5 over a range of b are shown in the same table. 
For the diffusive dimer model, we have formed [16,15], [15,16], [15,15], [16,14], [14,16], [15,14], and [14,15] Pade 
approximants to the z-series. The corresponding graphs are displayed in Fig. |9. It is seen that for small values of 
7, we obtain small estimates of 6, while for large 7, 8 — ► 1, suggesting the approach to the limiting saturation is 
via a mean-field like result, i.e. the i _1 power law. Hence we see that even though the exponential transformation 
Eq. ( |l4| ) works well for the exact series of of diffusive monomer model when 7 = 1/2, its use for general 7 is not very 
appropriate. We have also tried the transformation z = 1 — (1 + bt)~ x / 2 to the series for both diffusive models but 
the convergence is rather poor. 

We have tried and used a third method of extracting the saturation exponent 6. If we assume that for large enough 
times t, the saturation coverage 9 assumes a power law 

1-0 oct - *, (15) 

then we expect a plot of dln(l — 0)/dln(i) versus t or log 10 (f) should give a plateau of constant — S values. By forming 
[14,13], [13,14], [13,13] Pade approximants to the d\n(l — 9)/d\n(t) of the 28-term series for the diffusive monomer 
model with 7 = 1/2, we observe from Fig. [lO] that the agreement between different Pade estimates and the exact 
solution is excellent for log 10 (t) up to around 0.9. For diffusive dimer problem, three Pade approximants of [16,15], 
[15,16], and [15,15] are formed. The plots of dln(l — 9)/d\n(t) versus log 10 (i) for the diffusive dimer and monomer 
models, shown in Fig. [ll]and Fig. |lj, respectively, are obtained by taking the average of the 3 different Pade estimates. 
The graphs end before the difference between at least a pair of Pade estimates is more that 0.001. The last estimates 
in Fig. [ll] and Fig. |l2| are taken as the estimates for 5 and they are listed in the last two columns of Table |. These 
estimates for S are plotted in the same graph for the _F-transformed analysis for comparisons (Fig. ^ . It is seen that 
our last method of extracting the saturation exponents appears to be better than the F-transform analysis since it 
yields almost about the same estimates for 5. It does not involve any transformation which is not known in advance 
that will yield consistent results [E3| . Looking at the ends of the curves in Fig. (fl] and Fig. [|, we are certain that the 
power law regime is still not reached since the S estimates do not seem to converge to a constant value, except the 
case when 7 = 1/2 for the diffusive monomer model. From this we know that our estimates for 6 do not describe the 
true power law approach at large times t. Such information cannot be found in the F-transform analysis. We note 
that our last method of analyzing the series is easy to use compared to the _F-transform analysis. 



IV. MONTE CARLO SIMULATIONS 



To study the short and large time behaviors of the coverage, we have performed extensive and exhaustive simulations 
for the diffusive dimer and monomer models. For both models, we take an initially empty linear lattice with N = 20000 
sites with periodic boundary conditions so that the finite-size effects can be ignored. In each Monte Carlo step, a pair 
of adjacent sites is chosen randomly. The type of attempted process is then decided: deposition with probability p, 
where < p < 1, or diffusion with probability (1 — p). In the case of the deposition attempts, if any one of the chosen 
sites is occupied, the deposition attempt is rejected (unsuccessful attempt), else the adsorption attempt is accepted. 
In the case of diffusion, yet another decision is made either to move right or left, with equal probability. If the selected 
decision is diffusion to the right, we check the selected pair of sites are occupied and its right nearest neighbor site 
is unoccupied, then the dimer is moved by one lattice constant to the right. The left-diffusion attempts are treated 
similarly. In contrast to the diffusive dimer model, the diffusive monomer model allows monomers to move by one 
lattice constant. 

We define one time unit interval (At ~ 1) to be during which a deposition attempt is performed for each lattice 
site. Thus for N-site lattice, one unit time corresponds to N deposition attempts, on average. The diffusion rate 7 
relative to the deposition rate, is then 7 = (1 —p)/2p. 
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Straightforward simulation procedure, as described above, encounters a serious drawback in which at late stages, 
most adsorption and diffusion attempts are rejected. In order to study the behavior of the system at large times, we 
have used an event-driven algorithm to speed up the dynamics of the simulations 20 2^]. Let q be the probability 
that we can make a successful move, then the probability that the first (i — l)th trials is unsuccessful, and the ?th 
trial is successful is 



p< = q (i - q y- 



1,2,3,. 



(16) 



If we restrict all trials to be coming from the successful ones, then two consecutive trials are in fact separated by a 
random variable i in Eq. (|l6|). This distribution can be generated by 



ln(l - q) 



1, 



(17) 



where £ is a uniformly distributed random number between and 1. In employing this method, we have to keep and 
update an active list of successful moves/attempts, where from its length we can evaluate q at any instance. 

Simulations are performed on a cluster of fast workstations. Our numerical results are obtained for 7 = 0.05, 0.10, 
0.20, . . ., and 6.40, for t up to 2 20 . Each data set is averaged over 500 runs, and the longest run take about 150 
CPU hours on a HP712/60. The coverage (fraction of occupied sites), 0(t), is plotted in Fig. 13 and Fig. |lj for the 
diffusive dimcr and monomer models, respectively. We have also performed the simulation at 7 = 1/2 for the diffusive 
monomer model in order to compare the simulation results with the exact results. It is seen that the agreement 
between them are so good that actually an overlapping line is observed in Fig. [l4|. For 7 = 0, we have the exact 
solution M 



0(t) = 



1 - exp(-2(l - exp(-t))) 



(18) 



For extremely fast diffusion case, i.e., 7 = 00, exact results have been obtained, where 



t(d) 



1 / 1 



4 1 1 

for the diffusive dimer model 



l)-Iln(l-0), 



and 



-1 



(19) 



(20) 



for diffusive monomer model. The approach of (9 C — 9) ~ t^ 1 at all times is obvious for these extremely fast diffusion 
models. We have included the lines of slope —1/2 to indicate the i -1 / 2 power law clearly. It is seen from Fig. |l3| and 
Fig. [l4| that for 7 > 3.20, the system takes a very long time (t « 10 4 ) before it can enter the final i -1 / 2 regime. This 
explain why we have difficulty in extracting the actual power law approach from the primitive expansion in time t. 

To further confirm that the saturation approach indeed follows a power law, we have used a scaling analysis. For 
large times t, let us assume that (1 — 6) has the following scaling form 



\-6={ 1 ty^G( 1 b /t), 



(21) 



where G is a scaling function and b is a constant to be determined. Eq. (21) requires that G(u) tends to a constant 
when u tends to pi| . This is required because for large times t, (1 — 9) ~ t -1 / 2 . Let further assume that for large u, 
G(u) ~ u z for some constant z. For extremely large 7, we have (1 — 9) ~ t^ 1 (see Eqs. ( |l9|) and (|20|)) hence 2 =1/2 
and 6=1. Writing Eq. ( f2l| ) as 

\-e = (i/ 7 )( 7 A) 1/2 G( 7 A) 

EE (1/7)^/7), 

we then make log-log plots of 7(1 — 9) versus t/^y for the diffusive dimer and monomer models as shown in Figs. [l5| 
and 16, respectively. It is seen that the data collapse into single curves at large times t. Our numerical data confirm 
not only the power law decay but also a scaling behavior for (1 — 9). 
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V. CONCLUSIONS 



By using an efficient algorithm based on the hierarchical rate equations, relatively long series are obtained for two 
models of 1-d RSA with diffusional relaxation. Analyses of series are performed, but it is seen that even though the 
series are long, we only manage to extract the behaviors of the systems up to intermediate times only. To study the 
power law of a system at large times t, we see that a series which exhibits a continuous crossover behavior in its short 
and intermediate times ought to be long enough so that various orders of Pade approximants can still converge in the 
power law regime. Using these long series, we find that the analysis of series based on the ratio method by Song and 
Poland (25) was not useful. Specifically, for the diffusive monomer series when 7=1/2 (which corresponds to the A 
+ A — » process in Song and Poland's work), we obtain the saturation exponent (where they have used the symbol 
v) 5 = 2.0, 1.2, 0.895, 0.729, 0.624, 0.551, 0.498, 0.458, 0.426, 0.401, 0.380, 0.363, 0.351, . . ., which is not seen to be 
converging towards the expected value of 1/2. 

We have also performed extensive computer simulations using an efficient event-driven algorithm, where it allows 
us to use and simulate a larger system to much larger times t than it was done previously on a supercomputer [Q. 
The t~ x / 2 power law approach of 6 to its saturation is confirmed numerically at large times t. 
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FIG. 1. The data structures used to represent the rate equations Eq. (|3j) and Eq. (^). The first field of the node associated 
with 'o' is the representation of the pattern 'o'. The second field points to its first four derivatives (i.e. -2, 6, -22, 94). The 
third field is the highest derivative, h obtained so far for P{o); in this case h is 4. The rate equation is represented in the fourth 
field. The rate equation for the configuration 'oo' involves four configurations, one of them is 'oo' itself. 

FIG. 2. The plot of the inverted series of time t versus coverage 9 for the 28-term series with 7 = 1/2 for the diffusive 
monomer model. 

FIG. 3. The results for the saturation coverage 9 C as a function of the rate of diffusion 7 obtained from the DLog Pade 
analysis of the inverted series of t = t(9) with [16,15], [15,16], [15,15], [16,14], [14,16] Pade approximants for the diffusive dimer 
model. The saturation coverage estimates are very close to 1. 

FIG. 4. The results for the saturation exponent, S as a function of 7 obtained from DLog Pade analysis of the inverted series 
t — t{9) for the diffusive dimer model. These estimates for 3 are deduced from the residues associated with the poles in Fig. [| 

FIG. 5. The results for the saturation coverage 9 C as a function of the rate of diffusion 7 obtained from the DLog Pade 
analysis of the inverted series of t — t(9) with [14,13], [13,14], [13,13], [14,12], [12,14] Pade approximants for the diffusive 
monomer model. 

FIG. 6. The results for the saturation exponent, 5 as a function of 7 obtained from DLog Pade analysis of the inverted series 
t = t(9) for the diffusive monomer model. These estimates for 5 are deduced from the residues associated with the poles in 
Fig. I 

FIG. 7. The distribution of zeros and poles in the vicinity of (1, 0) for the [14, 13] Pade approximant for the 28-term series 
with 7 = 1/2 for the diffusive monomer model. A circle (cross) denotes a pole (zero). 

FIG. 8. Pade approximant estimates for the exponent 5, as a function the transformation parameter b, derived from the 
F-transform analysis of the 28-term series with 7 = 1/2 for the monomer diffusive model. 

FIG. 9. The plot of the values in Table |j. The diamonds and the squares denote the estimates of S from the F-transform 
analysis for the diffusive dimer and monomer models, respectively. The crosses and the circles correspond to the estimates of S 
from the dln(l — 9)/dln(t) versus log 10 (t) type analysis, for the diffusive dimer and monomer models, respectively. The error 
bars are smaller than the symbols and hence they are not displayed. 

FIG. 10. The analysis based on the dln(l — 9)/d\n(t) versus log 10 (t) plot for [14,13], [13,14], and [13,13] Pade approximants. 
We see that the agreement between three Pade estimates and the exact results is excellent for log 10 (i) up to around 0.9. 

FIG. 11. The plot of dln(l — 9)/d\n(t) versus log 10 (t) for the diffusive dimer problem. All curves end before the difference 
between at least a pair of estimates from [16,15], [15,16], [15,15] Pade approximants is greater than 0.001. The ends of curves 
for 7 = 0.1, 0.2, . . ., 1.0, 1.5, 2.0, . . . , 5.0, are displayed in the downward direction. 

FIG. 12. The plot of d ln(l — 9)/d \n(t) versus \og w (t) for the diffusive monomer problem. All curves end before the difference 
between at least a pair of estimates from [14,13], [13,14], [13,13] Pade approximants is greater than 0.001. The ends of curves 
for 7 = 0.1, 0.2, . . ., 1.0, 1.5, 2.0, . . . , 5.0, are displayed in the downward direction. 

FIG. 13. Monte Carlo simulation results for the diffusive dimer model. The sequence of 7 for curves between the exact curves 
for 7 = and 7 = 00, in the downward direction, is 0.05, 0.10, 0.20, . . . , 6.40. The line of slope —1/2 shows the t -1 ' 2 approach 
at large times t. 
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FIG. 14. Monte Carlo simulation results for the diffusive monomer model. The sequence of 7 for curves between the exact 
curves for 7 = and 7 = 00, in the downward direction, is 0.05, 0.10, 0.20, 0.40, 0.50, 0.80, 1.60, 3.20 and 6.40. Notice that 
the simulation results for 7 = 1/2 and the exact results agree with each other extremely well that only one line is seen. The 
line of slope —1/2 shows the i -1 / 2 approach at large times t. 



FIG. 15. The scaling plot for the diffusive dimer problem. The sequence of 7, in the downward direction, is 6.40, 3.20, 1.60, 
. . ., 0.05. The line of slope —1/2 is included to indicate the power law clearly. 

FIG. 16. The scaling plot for the diffusive monomer problem. The sequence of 7, in the downward direction, is 6.40, 3.20, 
1.60, . . ., 0.05. The line of slope —1/2 is included to indicate the power law clearly. 

TABLE I. The F-transform analysis gives the second and fourth columns which show the estimates of 8 deduced from the 
crossing regions of the graphs of 5 versus the transformation parameter 6, taken in the range indicated in the third and fifth 
columns. The last two columns show the results obtained from the dln(l — 8)/d\n(t) versus log 10 (t) type analysis. 
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TABLE 1. 



dimcr 



monomer 



7 



dimcr monomer 



0.1 0.431(3) 

0.2 0.499(2) 

0.3 0.542(2) 

0.4 0.573(3) 

0.5 0.599(4) 

0.6 0.623(5) 

0.7 0.649(4) 

0.8 0.669(5) 

0.9 0.685(5) 

1.0 0.702(5) 

1.5 0.759(5) 

2.0 0.792(5) 

2.5 0.830(3) 

3.0 0.854(4) 

3.5 0.87(1) 

4.0 0.885(5) 

4.5 0.898(4) 

5.0 0.907(3) 



0.269(4) 0.65 

0.395(4) 0.30 

0.441(2) 0.35 
0.47788(4) 0.45 
0.5061(5) 0.45 

0.535(1) 

0.557(2) 

0.576(1) 
0.5915(6) 
0.6050(6) 
0.6534(3) 
0.6822(4) 
0.7041(4) 

0.718(2) 

0.730(3) 

0.742(4) 

0.746(4) 

0.752(4) 



0.28 
0.30 
0.35 
0.35 
0.36 
0.38 
0.44 
0.47 
0.48 
0.51 
0.57 
0.58 
0.73 
0.80 
0.81 
0.88 
0.90 
0.95 



0.33 
0.35 
0.40 
0.40 
0.41 
0.43 
0.49 
0.52 
0.53 
0.56 
0.62 
0.63 
0.78 
0.85 
0.86 
0.93 
0.95 
1.00 



0.75 
0.85 
0.90 
0.95 
1.00 
1.25 
0.95 
0.75 
0.80 
0.80 
0.75 
0.85 
0.85 



0.70 
0.35 
0.40 
0.50 
0.50 
0.80 
0.90 
0.95 
1.00 
1.05 
1.30 
1.00 
0.80 
0.85 
0.85 
0.80 
0.90 
0.90 



0.407(1 
0.510(1 
0.566(1 
0.603(1 
0.618(1 
0.662(1 
0.682(1 
0.716(1 
0.673(1 
0.684(1 
0.810(1 
0.830(1 
0.836(1 
0.859(1 
0.865(1 
0.875(1 
0.884(1 
0.906(1 
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APPENDIX A: THE SERIES EXPANSION COEFFICIENTS FOR DIFFUSIVE DIMER AND MONOMER 

MODELS 

For the diffusive dimcr model, the series expansion coefficients P(o)(") for n = 0, 1, 2, ... , 31 are 



P(o)W 


= 1, 


P(o)« 


= -2, 


P(of) 


= 6, 


P(o)( 3 ) 


= -22, 


P(o)( 4 ) 


= 94, 


P(o)( 5 ) 


= -454-87, 


p(oy 6 > 


= 2430+ 136 7 + 32 7 2 , 


P(o)( 7 ) 


= -14214 - 16487 - 544 7 2 - I6O7 3 , 


P(o)( 8 ) 


= 89918 + 177767 + 6720 7 2 + 2752 7 3 + 8967 4 


P(o)( g ) 


= -610182 - 1836407 - 754887 2 - 348487 3 - 




-5376 7 5 , 



P(o) (10) = 4412798+ 18768247 + 8295047 2 + 406464 7 3 + 2078407 4 
+98560 7 5 + 33792 7 6 , 

P(o) (11) = -33827974- 192905607 - 9193152 7 2 - 4688896 7 3 
-2538496 7 4 - 1364672 7 5 - 651008 7 6 - 2196487 7 , 

P(o) (12) = 273646526 + 201202624 7 + 104153760 7 2 + 55017632 7 3 
+30829536 7 4 + 17641472 7 5 + 9613280 7 6 
+4487168 7 7 + 1464320 7 8 , 

P(o) (13) = -2326980998 - 2140434856 7 - 1212903168 7 2 
-664001632 7 3 - 381899008 7 4 - 226930976 7 5 
-132747296 7 6 - 71354592 7 7 - 31943680 7 8 - 9957376 7 9 , 

P(o) (14) = 20732504062 +23292327080 7 + 14540797344 7 2 
+8274623072 7 3 + 4870602688 7 4 + 2973218272 7 5 
+1818085760 7 6 + 1059060704 7 7 + 550524320 7 8 
+233132032 7 9 + 68796416 7 10 , 

P(o) (15) = -192982729350 - 259697613072 7 - 179401122720 7 2 
-106542051936 7 3 - 64143045632 7 4 - 40012549568 7 5 
-25252420800 7 6 - 15495356512 7 7 - 8810633856 7 8 
-4369323616 7 9 - 1734705152 7 10 - 481574912 7 n , 
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P(o)( 16 ) = 1871953992254+ 29690988160167+ 2275429063808 7 2 

+ 1416343409184 7 3 + 8726451 13632 7 4 + 554654077056 7 5 
+358789726592 7 6 + 2285941 11904 7 7 + 1379526891207 8 
+75460620608 7 9 + 35391538304 7 10 + 131050045447 11 
+3408068608 7 12 , 

P(o) (17) = -18880288847750- 34819585889272 7 - 29629807780320 7 2 
-194142844252807 3 - 12257502997152 7 4 - 79253532409607 5 
-52321902912967 6 - 34334002494087 7 - 2164670388352 7 8 
-12651632256327 9 - 6590408252167 10 - 2908435558407 11 
-1001936322567 12 - 243433472007 13 , 

P(o) (18) = 197601208474238+ 4188579227402167 +3956022991736967 2 
+273978158172160 7 3 + 177588283397088 7 4 
+ 116718001713792 7 5 + 78414883481952 7 6 
+527265513900487 7 + 34420850096512 7 8 + 21138424854304 7 9 
+118308654893447 10 + 5828807502880 7 n + 2414315305984 7 12 
+773310775296 7 13 + 175272099840 7 14 , 

P(o) (19) = -2142184050841734- 51673341163372487 

-5409353469693312 7 2 - 3974358277418464 7 3 
-2650710992073376 7 4 - 1770760359697088 7 5 
-12082501237828167 6 - 8294984596191047 7 
-5575141706110087 8 - 3563424407812487 9 
-210643968436256 7 10 - 1119959017880967 11 
-519510906260807 12 - 201796833657607 13 
-60136484372487 14 - 12707227238407 15 , 

P(o) (20) = 24016181943732414+ 65354875319253600 7 + 756740156950717767 2 
+59170426635099072 7 3 + 40709461342679712 7 4 
+276545533753431687 5 + 191391651825736327 6 
+ 133794779487778567 7 + 9220546462388512 7 8 
+6096023066059456 7 9 + 3768570926043360 7 10 
+2126161230748448 7 n + 1067953378669088 7 12 
+465021052820384 7 13 + 1694389076145607 14 
+470472447754247 15 + 92688010444807 16 , 

P(o) (21) = -278028611833689478- 847070521919796296 7 

-1082140708741735168 7 2 - 902807166330596704 7 3 
-6424631128778783687 4 - 4442021480352959367 5 
-3115570829540717767 6 - 221310190035972736 7 7 
-1558919759261379207 s - 106130622734975296 7 9 
-681632129900881607 10 - 404010832103840967 11 
-216285498791485127 12 - 102246856969797767 13 
-4170519535181536 7 14 - 1426800665301088 7 15 
-369849931661312 7 16 - 679712076595207 17 , 
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P(o) (22 '> = 3319156078802044158+ 112457240956831982807 

+ 15806297066594859744 7 2 + 14097805354882121664 7 3 
+ 10405357836443942528 7 4 + 7331412297306645056 7 5 
+5209126487226619872 7 6 + 3753693513361866688 7 7 
+2695682775965810496 7 8 + 1883121553855520800 7 9 
+ 1250352187435821024 7 10 + 7731192322639122567 11 
+436759047669796704 7 12 + 221031829068739616 7 13 
+98090146080366016 7 14 + 37415367315785664 7 15 
+12034887222344992 7 16 + 2918797259309056 7 17 
+500840477491200 7 18 , 

P(o) (23) = -40811417293301014150- 1528498191432711868647 

-235669362697015641504 7 2 - 2250304342568664740487 3 
-172730720539355533536 7 4 - 124209020387171547840 7 5 
-89396791275859021472 7 6 - 65265468008645791072 7 7 
-476812262349848156487 8 - 340789144796208485447 9 
-23302022034195791264 7 10 - 149501545807123630087 11 
-8846189597178252768 7 12 - 47458488757382929287 13 
-2265231509811056512 7 14 - 9419517996081430407 15 
-335416253872860704 7 16 - 101595547064760928 7 17 
-231071126996910087 18 - 37062195334348807 19 , 

P(o) (24) = 516247012345341914942+ 21258333457028609261287 

+35847097363169161971207 2 + 36676358259378862652487 3 
+29353022509972700083207 4 + 21579116585480016721927 5 
+15736120473774685552967 6 + 11627804829705368101127 7 
+862635165194788176384 7 8 + 629254138123646209984 7 9 
+4416704947841194150407 10 + 2927736974233096298887 11 
+180386574048628422176 7 12 + 101774434388607688384 7 13 
+51747976608776753696 7 14 + 23266775394437326144 7 15 
+9052145344807326144 7 16 + 30025168121887886087 17 
+8578155319416312967 18 + 183397309348839424 7 19 
+275319165340876807 20 , 

P(o) (25) = -6711185258405244576646- 302381800027045963332087 

-555982005708617079796807 2 - 60976389483506749681568 7 3 
-510033126453228642681287 4 - 384047404947730629501767 5 
-283899641877705734750407 6 - 21217507587703139716128 7 7 
-159596396322288053042567 8 - 11856781604498555126144 7 9 
-85199235046654863482567 10 - 58148140529469820132167 11 
-37130604095853993214727 12 - 2188983105017137174464 7 13 
-11754786171295814905607 14 - 566022591719541854336 7 15 
-239600608251046218528 7 16 - 870939539470198434247 17 
-26826301792524728832 7 18 - 7241192306564087872 7 19 
-1458610850414723072 7 20 - 2052379232541081607 21 , 
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P(o)( 26 ) = 89574471680939133937534+4396639783040107613093367 

+878864969592056661663680 7 2 + 1033220218894344066221 152 7 3 
+905162593426518959688480 7 4 + 699466982856423681626592 7 5 
+524550622363351486026048 7 6 + 396321792642538988104928 7 7 
+301868194425498598076992 7 8 + 2279892421914718162496967 9 
+ 1673333715204555998674247 10 + 1172432694506111524160007 11 
+772958577807145190609287 12 + 47370056124944928474784 7 13 
+266754532055173825258567 14 + 136278341853090747547207 15 
+62154481 1341 18540237767 16 + 2476949102568415569728 7 17 
+839899058850812411936 7 18 + 239165277641087918144 7 19 
+61092781558844261824 7 20 + 11620352278518562816 7 21 
+1534822730422026240 7 22 , 

P(o) (27) = -1226366187219563392423046- 65313799366796085554900487 

-141531585602297488311237127 2 - 178299425163416188756834887 3 
-163903886577137773843372487 4 - 13024188161033852626760832 7 5 
-9918027745419496427963072 7 6 - 7573935646469621224533504 7 7 
-5835396744531198025370240 7 8 - 4473355152862831429807584 7 9 
-3346869880800982736546432 7 10 - 24014960332843931063943047 11 
-1629516516241948614285760 7 12 - 10338019403116807545461127 13 
-607038025256747982659264 7 14 - 3264521288444553254076487 15 
-1587574024094091985719687 16 - 68639251401951128448448 7 17 
-257597449402422704669767 18 - 81333387647344437000647 19 
-21273771284716313801607 20 - 515040012051367812096 7 21 
-92704039977088974848 7 22 - 11511170478165196800 7 23 , 

P(o) (28) = 17208434165059531880467902 + 990815064203137109867727367 

+2321039034292190815329078727 2 + 313143052910337271943545632 7 3 
+3025318821374525912550772487 4 + 2476952038720645244693299207 5 
+ 191751044433982602086189184 7 6 + 148003121016823516846598336 7 7 
+ 115246426164717743454533376 7 8 + 89548623260586320065607520 7 9 
+68178700940595222195124032 7 10 + 499933834900064173166241287 11 
+34821933734970293328723264 7 12 + 22791830094404699697241952 7 13 
+13891323265031765512465344 7 14 + 78146658242265328529849607 15 
+4016982900101880080201440 7 16 + 1862205866467120855344800 7 17 
+764317776739694865201344 7 18 + 270291929328683129059200 7 19 
+792966295169514597200327 20 + 188795468379216957855047 21 
+43381109675360770634247 22 + 740404461632374177792 7 23 
+86564001995802279936 7 24 , 

P(o) (29) = -247289888972538586949878150- 15341726922400414526266365207 

-3874787071735125029536527360 7 2 - 5593914936389904746992944000 7 3 
-5687097195642308827551818432 7 4 -4806905879064545608139517504 7 5 
-3787785163974288880030601664 7 6 - 2955603108243908695631102784 7 7 
-2324464815975057231796395232 7 8 - 1828539431234039817220631104 7 9 
-1414541083036286620968100704 7 10 - 1058040639042272841296197632 7 n 
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-754795598374307863822508064 7 12 - 508234645790343970597871232 7 13 
-32032324581 14798314355688007 14 - 1875528278798677286825742407 15 
-1011940158947797542589169287 16 - 498087522562438325931462407 17 
-220571287898032207765959687 18 - 86105534305813552248551687 19 
-28719920025921023969999047 20 - 781079045618282247374272 7 21 
- 16717326284702027661 IOO87 22 - 365029516009751558996487 23 
-59189145767564206407687 24 - 6525593996606633410567 25 , 

P(o) (30) = 3636599975026505414628377086 + 242351928344885925550633814807 

+65825866944987465106280540640 7 2 + 1015881783229206313793466933767 3 
+1087925839779609198231643463047 4 +951063474381223214631612865287 5 
+76389237086435129240476170944 7 6 + 60282902962582074896568152544 7 7 
+47861156428195317146948964544 7 8 + 38077374671506940952543821152 7 9 
+298890648039810706025634258247 10 + 22767753081638239927853052832 7 n 
+166030997716152474500366471047 12 + 114731204180731476489940921607 13 
+7454372028423481067708753792 7 14 + 45239407843617679104900042247 15 
+2547674890838510347194642304 7 16 + 1321026477071424198196864544 7 17 
+6241589625370142481041382727 18 + 2647178897000485197608957127 19 
+98517274209783711927693664 7 20 + 31032886905928789132563200 7 21 
+7806681 142471034370768032 7 22 + 1477165355656235095757792 7 23 
+3068299363814525863769607 24 + 47352796340581207900160 7 25 
+4930448797436123021312 7 26 , 

P(o) (31) = -54690132113431117456486546054- 3904023063589740475544079980327 

- 1 137582424553489289447807997792 7 2 - 1874665509618062031 1152412912967 3 
-2116316617323806096783903601344 7 4 - 1916826452624015867653776105344 7 5 
-15716190630161851462586681341447 6 - 1255045527053951534543760304384 7 7 
-10056067174132445219605451348167 8 -8084112349675332189620711941767 9 
-6431160536491413353073062613447 10 -4981898222613772520297003313287 n 
-3707377011241920409491305626567 12 - 262369454324949422539656437184 7 13 
-1752686845167761477555475300807 14 - 109871960830604539165436865664 7 15 
-64284643706828638935486335392 7 16 -34893515783410832792353344544 7 17 
-17434380523486099607832839808 7 18 - 7930481323762069309980955744 7 19 
-3230836698378503824136216160 7 20 - 1149328018944519338194160160 7 21 
-342504050878779281408793504 7 22 - 79568648016407862847405792 7 23 
-13027718319821439160157152 7 24 - 2576298242354693706469632 7 25 
-379070964781663249235968 7 26 - 37330540894873502875648 7 27 . 

For the diffusive monomer model, we have 
P(o) (0) = 1, 

P(o) (1) = -2, 
P(o) (2) = 6, 
P(o) (3) = -22 + 47, 
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P(o)^ = 94- 4O7- I67 2 , 

P(o) (5) = -454 + 3167+ 1367 2 + 80 7 3 , 

P(o) (6) = 2430 - 2384 7 - 8407 2 - 6O87 3 - 4487 4 , 

P(o) (7) = -14214 + 181167 + 42407 2 + 30727 3 + 31367 4 + 26887 5 , 

P(o) (8) = 89918 - 1414327 - 139207 2 - 97287 3 - 131207 4 - 17664 7 5 - I68967 6 , 

P(o) (9) = -610182+ 1143564 7 - 520407 2 - 210727 3 + 164167 4 + 604807 5 
+1056007 6 + 109824 7 7 , 

P(o) (10) = 4412798 - 9606304 7 + 1860776 7 2 + 776864 7 3 

+366240 7 4 + 76544 7 5 - 285824 7 6 - 658944 7 7 - 732160 7 8 , 

P(o) (11) = -33827974 + 83906644 7 - 28877184 7 2 - 95896007 3 
-5139808 7 4 - 3368896 7 5 - 1464064 7 6 
+ 1317888 7 7 + 42465287 s + 4978688 7 9 , 

P(o) (12) = 273646526 - 761825992 7 + 377695696 7 2 + 89195520 7 3 
+45390432 7 4 + 35391872 7 5 + 28371968 7 6 
+ 15465472 7 7 - 5471232 7 8 - 280616967 9 - 343982087 10 , 

P(o) (13) = -2326980998+ 7184044444 7 -46546805367 2 
-6372498407 3 - 2881902087 4 - 2591729287 5 
-257515008 7 6 - 227319040 7 7 - 137749504 7 8 
+ 158412807 9 + 189190144 7 10 + 2407874567 11 , 

P(o) (14) = 20732504062- 702837112167 + 56240459224 7 2 
+2043271904 7 3 + 7315599367 4 + 12838916487 5 
+1667097280 7 6 + 1836321792 7 7 + 16842636807 8 
+ 1097334784 7 9 + 330752007 10 - 12965478407 11 
-1704034304 7 12 , 

P(o) (15) = -192982729350+ 7124956904687- 678540577872 7 2 
+439437292167 3 + 156897641607 4 - 248456832 7 5 
-6690721088 7 6 - 9906840832 7 7 - 112463311367 8 
-10574563328 7 9 - 76519096327 10 - 12065832967 11 
+9007038464 7 12 + 12171673600 7 13 , 

P(o) (16) = 1871953992254- 74749449904887 +82538992078087 2 
-1370040257024 7 3 - 3724111847687 4 
-943160904967 5 - 147027022727 6 + 15363239488 7 7 
+29522467968 7 8 + 35529220864 7 9 + 38287712768 7 10 
+40591720448 7 n + 15336308736 7 12 - 63292702720 7 13 
-87636049920 7 14 , 



32 



P(o) (17) = -18880288847750+810578141788607- 1017825061546647 2 
+27189637704816 7 3 + 53359581396807 4 + 1434188162432 7 5 
+658007059712 7 6 + 490637978560 7 7 + 4789881032967 8 
+504910140800 7 9 + 504004959232 7 10 + 347157901824 7 n 
-27525865472 7 12 - 1576231731207 13 + 4491347558407 14 
+635361361920 7 15 , 

P(o) (18) = 197601208474238- 907450604595520 7 

+ 1276490317628872 7 2 - 4698549342683207 3 
-57231645622432 7 4 - 125990101639047 5 - 81362015750407 6 
-88575623568007 7 - 10650704276864 7 8 - 12583232475008 7 9 
-142154519773447 10 - 146145044541447 11 
-114167769891847 12 - 3717889196032 7 13 + 14788583424007 14 
-32141810073607 15 - 46344005222407 16 , 

P(o) (19) = -2142184050841734+ 104759862861346447 

-163127694439152967 2 + 7646757249760992 7 3 
+364024229452512 7 4 + 209150481482887 5 

+609659980508167 6 + 104195753619392 7 7 + 1435293294746887 8 
+180347215997184 7 9 + 2140400194657287 10 
+2427792303614087 11 + 2514167247347207 12 
+201407890808064 7 13 + 80080105439232 7 14 
-132055498752007 15 + 231720026112007 16 
+339856038297607 17 , 

P(o) (20) = 24016181943732414- 124577072070506344 7 

+2126640858167030887 2 - 121030264397636800 7 3 
+3207001016515744 7 4 + 1926788775361664 7 5 
-97846534908224 7 6 - 10115252465879687 7 
-1613872175868224 7 8 - 2122232829566464 7 9 
-25824989142064007 10 - 3029609305192704 7 n 
-3468160142157312 7 12 - 36542161002588167 13 
-2984599157455872 7 14 - 1274373807013888 7 15 
+114396518154240 7 16 - 1681393031577607 17 
-250420238745600 7 18 , 

P(o) (21) = -278028611833689478+ 15244229655516791647 

-2830011720336190520 7 2 + 1893620809129623184 7 3 
-1903129004544258567 4 - 524295730139816327 5 
-5696486069346432 7 6 + 8870271163894976 7 7 
+165082584490355207 8 + 22318365235658944 7 9 
+273305067750203527 10 + 321275953226076807 11 
+37750764697213632 7 12 + 44585314916972224 7 13 



+487104587969351687 14 + 40826850701992704 7 15 
+ 18137503276204032 7 16 - 9712727831347207 17 
+1227059169853440 7 18 + 1853109766717440 7 19 , 
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P(o)( 22 ) = 3319156078802044158- 191767473598799236007 

+38453690823403414456 7 2 - 295616584530263056967 3 
+5155032487045587520 7 4 + 942243974334873216 7 5 
+ 1054392992597687687 6 - 811112039502624007 7 
-1609973662870471687 s - 215954557589465664 7 9 
-259399767952654784 7 10 - 2969561133721559687 11 
-3415104197362242567 12 - 415727829470043072 7 13 
-528023625074855232 7 14 - 615174223548672512 7 15 
-5349706099205982727 16 - 244776032213663744 7 17 
+8131502895267840 7 18 - 9000818866913280 7 19 
-137659582670438407 20 , 

P(o) (23) = -40811417293301014150+2477698929981489295407 

-5335512188401238714087 2 + 4631332184273136615687 3 
-1 1490803163561 1603008 7 4 - 130266516009884435207 5 
-6948357946361189767 6 + 985945653624581504 7 7 
+1559657950629139840 7 8 + 1927443716610738944 7 9 
+2162116688501650432 7 10 + 22864532476149870087 11 
+2406004053525926464 7 12 + 27933441718319714567 13 
+38678393652441 15904 7 14 + 57687570003084094727 15 
+74794587259394695687 16 + 68424471091129482247 17 
+3214950387582763008 7 18 - 67385809698816000 7 19 
+663268898321203207 20 + 102618961627054080 7 21 , 

P(o) (24) = 516247012345341914942- 32851190258773721809207 

+75591085206077738434887 2 - 7308871611544450832064 7 3 
+23509817276694078437767 4 + 1189059932972290110087 5 
-16790780767539633792 7 6 - 166134406726798796807 7 
-15548046461070182784 7 8 - 153849865199887199367 9 
-141614726489952351367 10 - 112763617162221640967 11 
-73240301726727911687 12 - 45874607262189486087 13 
-8052405463773808896 7 14 - 24820684069180820672 7 15 
-56960923428566366080 7 16 - 88298755351002966272 7 17 
-86378467061276064256 7 18 - 41668128760418795520 7 19 
+5541750393274368007 20 - 490786338216345600 7 21 
-7674113652110131207 22 , 

P(o) (25) = -6711185258405244576646+446610593131469882903807 

-109331606843615382157384 7 2 + 1164802789446980682623207 3 
-460117452703595525162567 4 + 4408507986046780216327 5 
+792105834451739504064 7 6 + 305738012971578476416 7 7 
+159590366941358300288 7 8 + 95740210483655793728 7 9 
+30567600232164220480 7 10 - 56955396275740956800 7 n 
-166253571659143480512 7 12 - 2815181305218803006087 13 
-3626134237014992529287 14 - 3265619828846365100167 15 
-57663159522038980416 7 16 + 4745022239191978737927 17 
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+ 1016510392155756708096 7 i8 + 1083954865492145854976 7 iy 
+537244083313897897984 7 20 - 4530850240533626880 7 21 
+3645203984752312320 7 22 + 57555852390825984007 23 , 

P(o) (26) = 89574471680939133937534- 622083509256483761778144 7 

+ 16139728320194499383274007 2 - 1877838552918714034411104 7 3 
+8791369069241047014024007 4 - 602031450836693343964807 5 
-204981100289787007365767 6 - 5129659282973757956032 7 7 
-15454272820002453521927 s - 81470942122801407680 7 9 
+ 12605551052897408107527 10 + 2895967959719246612096 7 n 
+48628707806632815766407 12 + 7053872326919183964416 7 13 
+9217493906833607526464 7 14 + 107504087107971329587207 15 
+103292461891841066681607 16 + 60798394022360304523527 17 
-24731006441416007582727 18 - 114239779271065282124807 19 
-135870102930179893882887 20 - 69251152065902548090887 21 
+368741160983891804167 22 - 271663623284698644487 23 
-43282000997901139968 7 24 , 

P(o) (27) = -1226366187219563392423046 + 88714630759920067384469967 

-243106648999963260408409607 2 + 306607937241291044016829767 3 
-165902933460410709169792327 4 + 21017359376278303366990087 5 
+414485759213945048620480 7 6 + 699970756874769711097607 7 
+97967196804432111909767 8 - 12989214479805658674880 7 9 
-342552242369675027628807 10 - 596778814332438967472647 11 
-895673168985101138005767 12 - 122844647469002814372800 7 13 
-158096743387017593763520 7 14 - 192783430108635063948928 7 15 
-2186989940598626647191687 16 - 2148089199293596740244487 17 
-150396704456781735005888 7 18 - 17357306803247632596736 7 19 
+ 125070331075601599340800 7 20 + 170681832575907668260864 7 21 
+895156623199069630627847 22 - 299007080543601819648 7 23 
+203092466220920733696 7 24 + 3262796998303316705287 25 . 
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